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Abstract 

We compute 2-loop QCD corrections to the hard coefficient functions which arise in 
the factorization formula for B — > decays in the shape-function region. Our 

calculation provides the last missing piece required for a NNLO analysis of inclusive 
semileptonic B decays, which may significantly reduce the theoretical uncertainty 
in the extraction of the CKM matrix element \Vub\- Among the technical aspects, 
we find that the 2-loop hard coefficient functions are free of infrared singularities 
as predicted by the factorization framework. We perform a brief numerical analysis 
of the NNLO corrections and include a discussion on charm mass effects. 
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1 Introduction 



Charmless semileptonic B meson decays are mediated by a 6 — -u quark transition. The 
study of inclusive and exclusive semileptonic B decays provide two independent paths 
for a determination of the strength of the flavour- changing interaction. In the Standard 
Model the quark transition is governed by the CKM matrix element Vub, which constitutes 
an important input parameter for many analyses in flavour physics. 

Since the determination of CKM matrix elements from exclusive semileptonic B meson 
decays requires the knowledge of non-perturbative form factors, the determination from 
inclusive semileptonic decays is a priori cleaner from the theoretical point of view. The 
total decay rate can be described by a local Operator Product Expansion (OPE) in a 
simultaneous expansion in Aqco/'^b and as{mb) [T]. The local OPE can be applied 
for inclusive B —>■ Xclv decays, where the non-perturbative information is constrained 
to some numbers, HQET parameters, which can be extracted simultaneously with \Vch\ 
from an analysis of the B meson decay spectra. 

A major complication arises in charmless semileptonic B decays where experimental 
measurements of the B Xuiu decay rate have to introduce kinematical cuts to suppress 
the B Xciiy background. This restricts the experimentally accessible information to 
the shape-function region in which the hadronic final state has large energy Ex ~ fTT't 
but moderate invariant mass Px ~ ""^feAqcD- The theoretical description of partial decay 
rates in the shape-function region is more complicated and gives rise to a multi-scale 
OPE in terms of non-local light-cone operators. 

The development of QCD Factorization [2] and Soft-CoUinear Effective Theory(SCET) 
[3] has widely improved our understanding of strong interaction effects in charmless B 
meson decays. A factorization theorem for the partial B —>■ X^C-v decay rate in the 
shape-function region has originally been proposed in the diagrammatic approach [4J and 
has later been formulated in the operator formalism within SCET. The factorization 
formula can be expressed in terms of the hadronic tensor W^^ , which encodes the strong 
interaction effects m B ^ X^ii^ decays, schematically as 



with hard coefficient functions Hij, a jet function J and a shape function S which de- 
scribes the distribution of the residual light-cone momentum of the 6-quark within the 
B meson [5] (the c's are some tensor coefficients). Whereas the hard functions and the 
jet function can be computed in perturbation theory, since they describe fluctuations 
with virtualities ~ ml and ~ Aqcn^b respectively, the shape function represents a non- 
perturbative input to the factorization formula. 

In this work we address perturbative corrections to the factorization formula ([1]). 
Next-to- leading order (NLO) corrections to the hard functions Hij and the jet function J 
have been computed in [6] and next-to-next-to-leading order (NNLO) corrections to the 
jet function have been worked out in The only missing piece for a NNLO analysis of 




(1) 
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inclusive semileptonic B decays in the shape-function region consists in the corrections 
to the hard functions which we consider in this work. 

A first step towards the computation of NNLO corrections to the hard functions has 
been taken recently in [8]. The authors of [8j computed 2- loop QCD corrections to the 
form factors which parameterize the h ^ u quark transition. We present an independent 
calculation of these 2-loop form factors and extend the analysis of [8J in two respects: 
First, we perform the subsequent matching calculation within SCET to extract the hard 
functions Hij from the form factors which are formally infrared divergent. Second, we 
keep a non-vanishing charm quark mass in our 2-loop calculation which has been set to 
zero in [S]. 

The 2-loop calculation that we present in this work is closely related to the calculation 
of NNLO vertex corrections in hadronic B decays which we considered in [9]. It is in fact 
somewhat simpler than the one in [9] since it involves only a small subset of Feynman 
diagrams and 2-loop integrals. A further complication arises in the analysis of [9] through 
the appearance of evanescent four-quark operators which makes the matching calculation 
rather involved. 

The outline of this paper is as follows: To set up our notation, we briefly recapitulate 
the necessary ingredients for the matching calculation in Section [2l The 2-loop calculation 
of the QCD form factors is described in Section [31 We perform the matching calculation 
in Section m where we present our results for the NNLO hard coefficient functions Hij. 
We extend our analysis to include charm mass effects in Section |5l We briefly analyze the 
numerical impact of the NNLO corrections in Section [6l before we conclude in Section [71 
We collect some results for the NLO and NNLO coefficient functions and for the charm 
mass dependent Master Integrals in Appendices [^IIUI 



2 Preliminaries 

The basic quantity for the analysis of strong interaction effects to inclusive semileptonic 
B meson decays is the hadronic tensor W^'^ , which is defined via the discontinuity of 
the forward matrix element of the correlation function of two weak interaction currents 
J/^ = uY{l - 75)&, 

W^""" = ^ Im {B{pB)\i j d^x e'«" T [f>'{0)r{x)] \B{pb)). (2) 

For the present analysis of inclusive semileptonic decays in the shape-function region it 
will be convenient to introduce light-cone coordinates. We introduce two light-like vectors 
n± with = and n^n_ = 2. Any four-vector can be decomposed according to its 
projections onto these light-cone directions and a two-dimensional transverse plane as 

= +P± + (n-p)^ = {n+p,p^,n_p). (3) 

In the shape-function region the momentum of the hadronic final state Xu scales as a 
hard-coUinear momentum, i.e. px ~ ^fe(l,A,A^) with = Aqco/^b- A factorization 
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theorem for the hadronic tensor W^'^ can be estabhshed in a two-step matching procedure 
QCD — i> SCET — i> HQET. In the first step hard fluctuations with virtuahty ~ ml are 
integrated out. The semileptonic current becomes to leading power in A 

with the B meson velocity v and x_ = {n^x)n_/2. The static 6-quark field /i^, is defined 
in HQET and the hard-coUinear w-quark field ^ and the Wilson line Whc in SCET. The 
momentum space Wilson coefficients can be obtained via 

Q{n+p) = j dse''^+'PQ{s). (5) 

In the second matching step hard-coUinear fluctuations with virtuality ~ Aqco^b are 
integrated out. This gives rise to another perturbative coefficient function, the jet func- 
tion J, and a remnant HQET matrix element, the shape function S. The factorization 
formula for the hadronic tensor becomes to leading power in A 

W^^ ^ ^ H,,{n^p) Tr {^'.{^^-^^ j duj J{pl) S{uj). (6) 

The hard functions Hij, which we consider in this work, are related to the Wilson coeffi- 
cients Ci of the semileptonic current by H^j = CiCj. We may therefore mainly focus on 
the matching relation (jl]) instead of the factorization formula for the hadronic tensor ([6]). 

The matching calculation simplifies if we work with on- shell quarks and regularize 
ultraviolet (UV) and infrared (IR) singularities in Dimensional Regularization (DR) (we 
write d = 4 ~ 26 and use an anticommuting 75 according to the NDR scheme). In 
this case the SCET diagrams are scaleless and vanish and the main task consists in the 
computation of the 2-loop QCD corrections to the left-hand side of (jl]). We parameterize 
the QCD result in terms of three form factors. 



{uip)\J'^\b{pb)) = u{p) 



F,{q')r{l - 75) + F,{q')^{l + 75) + F,{q')^{l + 75) 

nib 



KPb), 
(7) 



with momentum transfer q = pb — p. On the other hand we work in a frame with 
V = [ri- + n+)/2 and choose the basis of three independent Dirac structures in (jl]) as 

15^ = 7^(1-75), r^ = t;'^(i + 75), r^ = n(i(i + 75). (8) 

Writing the momentum of the on-shell 6-quark as pb = mbV and the one of the w-quark 
as p = {n^p)n^/2 provides the link between the form factors Fi and the coefficient 
functions Ci. We finally introduce a dimensionless variable u to parameterize the light- 
cone momentum of the M-quark as m = n^p/irib. In this notation the momentum transfer 
becomes q^ = uml with u = 1 — u. 
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3 NNLO calculation of QCD form factors 



3.1 2-loop calculation 

To NNLO the QCD calculation of the form factors Fi{q^) involves the calculation of the 
2-loop diagrams from Figure [H As pointed out in the introduction, the calculation can 
be inferred from the calculation of the 2-loop vertex corrections to hadronic B decays 
The present calculation involves only a small subset of 2-loop diagrams and about half 
of the Master Integrals (Mis) that have been computed in [9]. 

The calculation follows the same strategy that has been used in [9]. We first express 
tensor integrals in terms of scalar integrals with the help of a general tensor decompo- 
sition. We then use an automatized reduction algorithm, which is based on integration- 
by-parts techniques [TU] and the Laporta algorithm [TT], to express the scalar integrals 
in terms of an irreducible set of Master Integrals (Mis). As long as we neglect the charm 
quark mass, the present calculation gives rise to 14 Mis which are depicted in Figure [2l 
All Mis have already been calculated in [9] with the help of the method of differential 
equations [12], the formalism of Harmonic Polylogarithms (HPLs) [^13], Mellin-Barnes 
techniques [H] and the method of sector decomposition [T5]. 

The Mis from Figure [2] can be expressed in terms of HPLs of weight w < 4. Through- 
out this paper, we give our results in terms of the following minimal set of HPLs 

i/(0;x) = ln(a;), i7(0, 0, 1; a;) = Li3(a;), 

if (1; x) = - ln(l - x), if(0, 1, 1; x) = 81,2(0;), 

H{-1; x) = ln(l + x), H{0, 0, 0, 1; a;) = Li4(a;), 

H{0, 1; x) = Li2(x), H{0, 0, 1, 1; x) = 82,2(0;), 

H{0,-l;x) = -U2{-x), /f(0,l,l,l;x) = 8i,3(x). (9) 

Moreover, we introduce a shorthand notation for a HPL of weight w = 3, whose expression 




Figure 1: 2-loop diagrams. Diagrams that result from mirroring at the vertical axis 
are not shown. The bubble in the last diagram represents the 1-loop gluon self-energy. 
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Figure 2: Scalar 2-loop Master Integrals. We use dashed (double) lines for massless 
(massive) propagators. Dashed/ solid/ double external lines correspond to virtualities 
O/nm^/m^, respectively. Dotted propagators are taken to be squared. 

in terms of Nielsen Polylogarithms is rather involved 

Hiix) = H{-l,0,l;x) 

= Lis (i^) + Li3 ([^) + S^M + Si,2(-x) - Si,2 (^) 
+ ln(l - x)U,{x) - In (^)Li2(-x) - In (^)Li2(^) 

and for a HPL of weight w = 4, which cannot be expressed in terms of Nielsen Poly- 
logarithms and has to be evaluated numerically, 

7^2(2:) =i/(0,- 1,0,1; x) = r^x'^^i^. (11) 

Jo ^ 

The Mis from Figure [2] have been computed in an analytic form [9] , apart from a single 
number Cq which arises in the calculation of the boundary condition to the 6-topology MI 
(last diagram from Figure[2]). With the help of the Mathematica package FIESTA [16j, 
we can compute this number to higher precision than in We now find Co = —60. 250 ± 
0.001. 

3.2 Renormalization 

The form factors, which we obtain from the 2-loop diagrams of Figure [T], are UV- and 
IR-divergent. The UV-divergences are subtracted in the renormalization procedure which 
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amounts to the calculation of several 1-loop diagrams with standard QCD counterterm 
insertions. We account for the wave-function renormalization of the quark fields by 

where Z2^u (-^2,fe) is the wave-function renormalization factor of the massless w-quark 
(massive 6-quark). Writing the perturbative expansion in terms of the renormalized 
coupling constant 



SI 

oo 



k k 



k=0 k=l 

we find that the renormalized form factors are given to NNLO by 

i j,bare' 

_ pW , 1 / 7(1) , 7(l)^ f(0) 

-^i ~ -^j.bare ~'~ 2 I ) "^i.bare' 

p{2) _ p(2) 1 / 7(1) p(l) /I / 7(2) , 7(2)N 1/7(1) yW\A p 

- ^i,bare + 2 I ^2,« + ^2,6 I ^i,bare + I 2 ^.^ + ^2,b ) " g 1^2,m " ^2,6 j I 



bare ' 



(14) 

We renormalize the coupling constant in the MS-scheme (adopting the notation f|T3|) ) 



z!,'' = -{^CA-"-n,T,]^_, (15) 



11^ 2 \1 

6 3 ^ J e' 

whereas the quark wave-functions and the 6-quark mass are renormalized in the on-shell 
scheme, 

and Z2^l = because of the fact that scaleless integrals vanish in DR. According to ffT^ . 

we also need the 2-loop expressions for the wave-function renormalization factors. For a 

(2) 

massless quark Z2 „ can be calculated easily, since it receives a contribution from a single 
diagram, shown in Figure [3l which introduces a mass scale. We find, in agreement with 

e.g. m, 



^ ml J (1 -£)(2-£)(l + 2£)(3 + 2£) ^ ' 




Figure 3: 2-loop diagram with massive fermion loop. 
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For a massive quark zfl has been calculated in [T7] . 



7(2) _ r 



+ 



'9 11 

^Cf + -Ca - 2niTF 



4 



12 



11 

y 



+ Cf - + ^riiTF + ( 1 + 4l]Tf 



+ " ^ + 16vrMn(2) - 24C3 + |l + QL^ ) 

+ ( + Svr^ - Svr^ ln(2) + I2C3 -^^L- ^L'' ] C 

, 113 47r2 38 , .9 , ^ 
+ ( — + — + y + ) n^Tp + 



6 2 

18 3 I ^ K ) 



(18) 



with L = In fJ, /ml and the number of massless quarks rii = 4. 



3.3 Renormalized form factors 

We now present our results for the renormalized form factors Fi{q^). It should be noticed 
that these form factors are formally IR divergent. We address the subtraction of these 
IR-divergences in the following section. 

It will be convenient to further decompose the coefficients F^^^ of the perturbative 
expansion (fT3|) according to 

^ = E (19) 
3 

We suppress an overall pref actor igwVub/'2\/^ and recall that = uml with u = 1 — u 
and L = In/x^/m^. In our normalization the form factors become in leading order 

Fi'\u) = l, FfH = 0, Ffiu)=0. (20) 

In NLO we give our result up to terms of 0{e'^), which are required for the IR-subtractions 
that we consider in the following section. In terms of a set of coefficient functions fi{u), 
2 = 1, .., 10, which we list in Appendix [XI we find that the form factor Fi is given by 



Ff%{u) = ~Cf, 

F}%{u) = CF(^fi{u)-Ly 



Fi^iu) = CF{f2{u) + Mu)L-h' 
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F«H = (^Mu) + Mu)L + ^Mu)L' - h'^ , 

FgH = (^U{u) + h{u)L + \f2{u)L' + \h{u)L' - ^L^) , 
whereas the other two form factors are finite at NLO and read 

F^]l{u)^CF[f,{u) + U{u)L^, 

fQ{u) = Cf(^Mu) + h{u)L + \h{u)L'^ , 

F^^l{u)=Cp(^h{u)+h{u)L^, 



(21) 



(22) 



and 



(23) 



In NNLO the divergent terms of the form factors can also be expressed in terms of the 
1-loop coefficient functions fi{u). The divergent terms of the form factor Fi read 



Fi%{u) = C'f(^L- f,{u)^ + ^CaCp - uiTpCf, 



L^-2h{u)L + -h{uf-h{u) 



+ CaCf 
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L-h{u) 



67 TT^ 

36 ^ 12 



F?Mu)^Cl 
-6C3 



+ uiTfCf 
2 



+ -LTfCf, 



-L' - 2h{u)L' - ( 2/2(1.) - h{uY ]L + h{u)h{u) - h{u) 



3 TT^ 



+ CaCf 



67\ / ^ ^ \ 461 ITtt^ 11/ 



+ uiTfCf 



10 



9 



-{L-n{u)\-- + - 



25 TT^ 



54 6 



TT 



2L'--Mu)L + - 



(24) 



whereas the other two form factors start with a 1 /e^-singularity, 



7(2) 
2,-2 



u 



-Clh{u) 



F?l{u) = CU h{u)h{u) - U{u) - 2U{u)L 



(25) 



and 



-Clh{u) 



(26) 



The finite terms of the 2-loop form factors give rise to a new set of coefficient functions 
ki{u), i = I, .., 10, which we collect in Appendix [Bl We find 



-L' - -h{u)L' - 2Mu) - h{uY ]L'-[ 2h{u) - 2h{u)h{u) + - 



+ CaCf 



/ 11 . / N /t)Y 7r-\ ^ , , 4t)i iYvr- , -, . \ . , . 



+ niTpCp 



'67 



— ^ + —fiiu) \ 

18 V 6 ^ ^ 18 6 ' 



461 177r^ 



^L' - ( 9/1 («) - Y - ( i^/2(n) + ^/i(n) + _ _ )L 



20 



25 _7r2 
27 ~ ^ 



TpCp 



F?i{u)=Cl 



14 /4 27r^ 

-L^^ - 2/i(«)L2 - i-Uu) ~—]L + Hu) 



+CaCf 



h{u) + 2 Mu)Uu) - U{u) \L- 2h{u)L 



h{u) + y/5(M)L 



+ niTpCp 



4 14 4 ■ 

3/6(m) - y/5(m) - ^h{u)L 



+ TfCf 



kiiu) - -h{u)L 



Fn{u)=Cl 



h{u) + 2 h{u)h{u) - h{u) ]L- 2U{u)L 



+CaCf 



+ niTpCp 



+ TfCf 



8 4 14 4 ■ 

— ln(M) - -h{u) - —h{u) - -h{u)L 



4 



kwiu) - -h{,u)L 



(27) 



We compared our results for the 2-loop form factors with the ones given in [8] and found 
complete analytical agreement. 
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4 Matching calculation 



4.1 IR-subtractions 



In order to extract the Wilson coefficients Cj from the form factors Fi of the previous 
section, we have to consider the perturbative expansion of the hadronic tensor W^'^ and 
to absorb the IR-singularities of the form factors into the jet function J and the shape 
function S. As the convolution of jet and shape function is independent from the indices 
(i, j) it is convenient to rewrite the factorization formula 



as W"" = Y.c'i-W^j with 

(2^ 



where the symbol ® represents the convolution integral. When we write the perturbative 
expansion in terms of the renormalized coupling constant, we should keep in mind that 
the actual 2-loop calculation has been performed in QCD with five active quark flavours, 
while the NNLO calculation of the jet function [7] and the partonic shape function [18] 
have been performed in SCET with four active flavours. We therefore have to account for 
an additional finite renormalization to consistently express the perturbative expansion in 
terms of the renormalized coupling constant of the four-flavour theory. 



(4), 



(29) 



with (cf. e.g. [19]) 



4, /2 . 7r2\ /2 . 4 



(30) 



In the four-flavour theory the perturbative expansion of the hadronic tensor becomes 

= (c«cf + [J ® s] + c^'^cf [J ® S] 

y^\f = (cf ^C'f + cf + Cl''>Cf - 6af^C\'^Cf - 6a^pcf^cf) [J ® S] 

+ (c\'^Cf + Cf^Cf^ - <5a«Cf ) [J ® S] + Cf ^Cf [J ® S] (31) 

As we have mentioned in Section [2l the matching calculation simplifies in our setup due 
to the fact that the SCET diagrams are scaleless and vanish. The IR-subtractions are 
therefore entirely determined by the counterterms. At NLO the counterterm can be 
inferred from [6]. In terms of the coefficient function fi{u) from Appendix [Al it reads 



[J®S] 



(1) 



L-h{u] 



[J®S] 



(0) 



(32) 
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In NNLO the counterterm can be extracted from the analysis in [TJ [T2]. We obtain 



[J®S] 



12C3 \Cf 



+ 



461 ITtt^ 



108 12 
TT^ 25 20 



4.2 Hard coefficient functions in NNLO 



[J®S] 



(0) 



(33) 



We now have assembled all pieces required for the NNLO calculation of the hard coeffi- 
cient functions Hij. Performing the IR-subtractions as described in the previous section, 
we observe that the UV- and IR- divergences cancel in the hard coefficient functions which 
represents both a non-trivial confirmation of the factorization formula ([6]) and a stringent 
cross-check of our calculation. 

We give the results for the hard coefficient functions in terms of the Wilson coefficients 
Cj, which arise in the matching relation (jl]) of the heavy-to- light current, Hij = CiCj. 
We remind that the perturbative expansion is formulated in terms of the renormalized 
coupling constant in a four-flavour theory and that in leading order 



0, 



0. 



(34) 



In NLO we obtain in terms of our results from Section [3l3l 

i^SM- (35) 
The NNLO coefficient functions were so far unknown. They are found in this work to be 



C?\u) = Fi'>+TFCF 



(2) 



^Cs + y/i(«) + ^(6/2(«) - 7r')L + 2h{u)L^ - 



U{u) - h{u)h{u) + (2h{u) - h{u)h{u)]L 



(36) 
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and 



F^F 



;h{u)L 



+ cl 



h{u) - h{u)h{u) + ( 2U{u) - h{u)h{u) ]L + -h{u)L^ 



(37) 



and 



-Cf\u) = F^i^+T,C, 



u 



\h{u)L 



+ cl 



Mu) - fiiu)Mu) + { 2Mu) - Mu)fs{u) )L + -h{u)L' 



(38) 



The hard functions Hij and hence the Wilson coefficients Ci are renormahzation scale 
dependent. According to the factorization formula ([6]), this scale dependence has to cancel 
against the one of the convolution of jet and shape function. The scale dependence of 
the Wilson coefficients Ci is governed by the renormahzation group equation, 



d 



dlnn 



Ci{u;fj,) 



rcusp(as)ln h7 [as) 



Ci{u;iJ,), 



(39) 



where Fcusp is the universal cusp anomalous dimension which describes the renormahza- 
tion properties of light-like Wilson lines [20]. Writing the perturbative expansion as 



Tcusp [^s) 



CO , CO , 



(40) 



k=l k=l 

we derive from our NLO expressions in ( 135]) the familiar results 
Similarly, we can read off from our NNLO results (|36l) - ( |38l) that 



(41) 



r(2L = CaCf 



cusp 



268 4tt^ 



9 3 
-^-24<3 



—niTpCF, 



+ CaCf 



22C3 



1549 77r2 



54 



6 



+ niTpCp 



250 27r2 



(42) 



While our result for the 2- loop cusp anomalous dimension is in agreement with [20], the 
expression for 7'*^^^ has not yet been computed directly so far. Our result con&ms the 
conjecture in [2JJ, which was based on known results of the 2-loop anomalous dimensions 
of the jet and the shape function. 
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5 Charm mass effects 



In this section we extend our analysis to include charm mass effects. The charm quark 
enters the 2-loop calculation through the 1-loop gluon self energy in the last diagram of 
Figure [H So far we treated the charm quark as massless and the respective effects are 
given in the above formulas by the terms that are proportional to rii = 4, which denotes 
the number of massless quarks. 

A comment is in order about the power counting of the charm quark mass within 
the effective theory approach. One possible choice is given by ~ jJhc ~ (AqcD^^b)^''^, 
i.e. the additional mass scale is considered to be of the order of the intermediate hard- 
collinear scale. In this case charm mass effects enter in the calculation of the jet function. 
In this setup the charm mass represents an IR-scale in the QCD— s>SCET matching cal- 
culation and can therefore be set to zero. If we adopt the scaling rric ~ fihc, the hard 
coefficient functions Hij are thus given by the formulas from Section HI 

A second choice, which has been widely applied in studies of 6 — > c decays, is to 
consider the charm quark in the heavy quark limit with rrib oo and rric ^ oo while 
mc/mb is kept fixed. In this scenario the hard coefficient functions Hij depend on the 
mass ratio z = rric/rnh. On the other hand the jet function does not dependent on the 
charm quark mass in this setup, since it is defined in SCET with three active quark 
flavours. Both approaches differ in the fact that the first one allows to resum logarithms 
of the form hmtc/mb while the second one does not. 

In the following we consider the second scenario, i.e. we keep a finite charm quark 
mass in the 2-loop QCD calculation and the subsequent matching calculation. We briefly 
outline the steps of the calculation that change when we introduce a finite charm quark 
mass and write our results schematically as (for an arbitrary quantity Q) 

Q{z) = Q + AQ{z), (43) 

where Q refers to the respective expression from Section [3] or H] but with ni = ?> massless 
quark flavours and AQ{z) gives the additional contribution from a massive charm quark. 

5.1 QCD form factors 

The charm mass enters the QCD calculation through the 2-loop diagram with a closed 
charm quark loop, which gives rise to four additional Mis (depicted in Figure Hj). Since 
these Mis have not been considered in [9j, we present our results for these Mis in Ap- 
pendix O As discussed in the appendix, we obtain analytical results apart from the finite 




Figure 4: Additional Master Integrals with propagators of mass rUc (wavy line). 
Notation as in Figure\^ 
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terms of the 4-topology Mis. We could, however, find suitable representations which al- 
low us to evaluate these terms numerically to very high precision. Moreover, we provide 
compact parameterizations for physical values of z = rric/mfy in the appendix, which 
reproduce our numerical results at the percent level. 

The charm mass also modifies the renormalization procedure. We get an additional 
contribution to the 2-loop wave-function renormalization factor of the u-quark. 



AZi%)=TFCF 



z'^ml 



Tie) 



2e(3-2e)(l + e) 



;i -£)(2-e)(l + 2£)(3 + 2e)^ 



whereas the one of the 6-quark receives the additional contribution [T7] . 



AZg (z) = TpCpi ( 1 + 4L J - + + ( ^ + 8 In(^) )L, + 



22 



(44) 



(5 - 18z - 30z^ + 12z^)7r2 



+ -(19 + 12^2) Ij^(^) ^ 8(^1 _^ 3^4^) j^2(^^^) _ g^(^3 ^ 5^2^) j U2{-z) + \n{z) ln(l + z] 



z){2- z- z^ -6^) Li2(z^) + 21n(z) ln(l - 2^) 



443 



28^2 + 0{e] 



(45) 

where we introduced Lc = In /i^/m^ = L — 21n(2;). 

These changes modify our NNLO expressions for the renormalized QCD form fac- 
tors Fi{q^) from Section 13. 3[ We now introduce three additional coefficient functions 
ki{u,z),i = 11,. .,13, which we list in Appendix [Bl The form factor Fi receives the 
additional contribution 



AFl%iu,z) 
AfI'Uu,z) 



AF[%{u,z) 



4 



-LcTfCf, 



2LI-^[Mu)-2Hz))l^ + ^ 
Ll-2(^Mu)~2\n{z))Ll 
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- - Uu) + 2 ln(z)/i(n) - 2 In^(^) -—]L, + k^iu, z 



6 



TfCf, (46) 



whereas the other two form factors are modified by 



AF^i{u,z) 



AF^^{u,z) 



ki2{u,z) - -f^{u)L^ 



TpCp- 



(47) 
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5.2 Hard coefficient functions 



According to the power counting that we adopt for the charm quark mass scale, the charm 
quark is integrated out in the QCD^SCET matching calculation. The IR-subtractions 
turn out to be analogous to those of Section 14.11 apart from the fact that we now match 
onto SCET with three active quark flavours. If we express the perturbative expansion in 
terms of the renormalized coupling constant of the three-flavour theory, we thus have to 
account for an additional contribution to (129!) with 



ASai'\z)=TF 



4r /2^2 

-Lc+ -Lt H ]e 

3 V3 " 9 



(45 



where we recall that Lc = In/x^/m^ = L — 2\n{z). 

We summarize our results for the hard coefficient functions Hij = CiCj. If we keep 
a finite charm quark mass in the matching calculation, which corresponds to the power 
counting rrib —>■ oo and rric —>■ oo with rric/mb fixed, we find that the coefficients functions 
are given in NNLO by our results from (136|1 - (l38l) with rii = 3 massless quarks and the 
additional contributions 



h,iu,z)+{Mu)~2\niz)]^ + ^C3 



ACi^\u,z) = ku{u,z)TFCF, 



ACP{u,z) = kn{u,z)TpCF, 



u 



(49) 



which depend on the quark mass ratio z = rric/mi,. 



6 Numerical analysis 

We briefiy discuss the numerical impact of the considered NNLO corrections. As a 
phenomenological analysis of inclusive semileptonic B decays is beyond the scope of 
the present paper, we illustrate our results at the level of the Wilson coefficients Ci 
which arise in the matching relation (jl]) of the semileptonic current (recall that the hard 
coefficient functions from are given by Hij = CiCj). We in particular do not study 
the renormalization scale dependence, since the Wilson coefficients are explicitly scale 
dependent, cf. ( 1391) . The scale dependence is cancelled only at the level of the hadronic 
tensor against the one of the jet and the shape function. 

In Figure [5] we show the Wilson coefficients Ci at the scale fi = mi,. We see that 
the considered NNLO corrections give moderate contributions to the Wilson coefficients, 
which add constructively to the NLO corrections. The effect of the finite charm quark 
mass is small, but amounts to ~ 10 — 20% of the NNLO contribution as can be seen in 
Figure [61 In Figure [6] we also illustrate a scenario, which considers the charm quark to 
be as heavy as the 6-quark, i.e. z = rric/mh = 1 (dashed line). Somewhat surprising, the 
curves for physical charm masses with 2; ~ 0.3 are much closer to the limit rric = rrib 
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0.01 r 




Figure 5: Wilson coefficients Ci at the scale n = nih as a function of u (the mo- 
mentum transfer is given by = (1 — u)m'^). The dashed lines refer to the NLO 
results and the solid lines to the NNLO results with z = mc/mh = (orange/light 
gray) and z = 0.3 (blue/dark gray). 
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Cf\u) 



-0.05 
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-0.15 
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0.010 r 



0.008 



0.006 - 



0.004 - 



0.002 




-0.005 



-0.010 



-0.015 



-0.020 - 




Figure 6: NNLO contributions to the Wilson coefficients Ci at the scale fi = as 
a function ofu. The solid lines again refer to the NNLO results with z = nic/mi, = 
(orange/light gray) and z = 0.3 (blue/dark gray). The dashed green/light gray line 
shows the limit z = 1. In the first and the third plot the dashed line can hardly be 
distinguished from the solid blue/dark gray line. 
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than to the hmit = (notice that the effect is even shghtly larger for physical charm 
masses than for rric = mf,). 

7 Conclusion 

We computed NNLO corrections to the hard coefficient functions Hij which arise in the 
factorization formula (jS]) for inclusive semileptonic B meson decays in the shape-function 
region. Together with the 2-loop corrections to the jet function from [7], our calculation 
completes the NNLO calculation oi B ^ Xud-v decays in the shape-function region. 

The considered calculation is closely related to the 2-loop calculation in charmless 
hadronic B decays [9J. In particular, all Master Integrals that appear in the present 
calculation have already been calculated in [9]. We confirmed analytical results for the 
2-loop form factors which have been computed recently in [8] and extended the analysis 
to extract the hard coefficient functions Hij from the form factors which are formally 
IR-divergent. We in addition kept a finite charm quark mass in the 2-loop calculation. 

A phenomenological analysis of partial decay rates in the shape-function region was 
beyond the scope of the present paper. We briefly discussed the numerical impact of the 
considered corrections to the short- distance coefficients, which was found to be moderate. 
One should keep in mind, however, that the present calculation reflects only a part of the 
full NNLO calculation. We in particular expect that the renormalization scale dependence 
of the theoretical prediction for partial decay rates will be significantly reduced, once the 
renormalization group improvement is taken into account in NNLO. 

Note added: While this paper was in preparation, the work in [221 [23] on the same 
topic appeared. While the analysis in [23] did not go beyond the one in [S], the authors 
of (22] also performed the full matching calculatioiil|. We have compared our results 
for the NNLO hard coefficient functions from fl36|) - fl38|) with those in [22] and found 
agreement. None of the works in [221 [23] addressed charm mass effects. 
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A NLO coefficient functions 

The 1-loop calculation of the QCD form factors gives rise to the following coefficient 
functions 

/2N = -6 - Ti^ ^ln(«) -2ln'{u) - 2U,{u), 

12 u 

, , . 57r^ 1 ^ 24(1 - 2(7) - tt^m , , , 1 - 3-u / , o/ x t • / x\ 
h{u) = -12 - — + -C3 ^ InH + ^—(\Ti\u) + Li2(M)j 

+ ^ln3(M) +41n(M)Li2(M) -2Li3(«) +4Si,2(m), 

vr^ 5 96(1 - 2iZ) + (1 - 3n)7r^ + SfZCa , , , 

fA\u) = —24 — Cs ^ mm 

^ ' 2 160 6^ 12u ^ ' 

^ 24{l-2u)-n'u ^^^, + Li2(^)) - ^n^(K) - 41n^(K)Li2(K) 

- ln3(ii) + 2 ln(^i)Li2(iZ) - Usiu) + 281,2(1*)^ - 8 \n{u)Si,2{u) 

+ 41n(ii)Li3(iZ) - 2LU{u) - 8Si,s{u) + 4S2,2{u), 

2 2u 
h{u) = - + — ln(M), 



U V? 



h{u) = - + ^(ln(..) - ln^(..) - Li2(«)), 

^ , , 48 + 7r2 2M/24 + 7r2 , , ^ ^. . ^. , , 2 , o , , r.^ r ^T ■ r ^ 

fiiu) = g_ + ^ ( ln(M) - ln2(M) - Li2(m) + - ln='(M) + 2 ln(M)Li2(M) 



Li3(M) + 2Si,2(M) , 



2 2(1 - 2^) , 



4 2(l-5«), , , 2(1 - 2m)/ 2/ X T•/-^ 
/g M = -- - — ^ In M + — \n\u) + L12 M 



^ , , 48 + 7r2 24(l-5M)+7r2(l-2ii), , , 2(1 - 5m) / , 2, x , t • /-n 
fio{u) = — ^ ^3^-^ ^ln(M) + ln2(M) + Li2(m) 



6m 6m^ u 



2(1 - 2m) /2 



M^ 



Q ln=^(M) + 2 ln(M)Li2(M) - LigCiZ) + 2Si,2(fZ)^ . (50) 
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B NNLO coefficient functions 



In the 2-loop expressions of the QCD form factors (!27l) we introduced the coefficient 
functions ki{u), which read 

2{5 + 31u + 3u^ + 3u^)^. 4(13 + 53m + 15^2 + 3m3)^ 
ki{u) = Li4(m) H S2,2(m) + 8Si,3(m) 



2(3 + 14m + 3^2 + #) , , X. , X , / / X 16 

^ ^ ^Li2 M ' - 8 In u Lis u + 161n m Si,2 m + — 

u'^ 3 



, ^r^ 2, ^, ■ i -^ , 3 + IOOm + 145^^ + 4m\ . , , , 42 - 68m , 3, , 
+ 161n'(M)Li2(M) + —- Li3(u) + — In^(M) 

2(3 - 242m- 233^2 -32m3)^ 2(6 + 41m + TTm^ + 2m3) , ^ . 
^^'^(^^ + "^^u l-WLi2(^) 

3 + 89m + 785m2 + 635m3 (5 + 31m + 3m2 + 3M3)7r2\ ^ ^ 

T^^^ + ^ 3 — Li2 n 

18m"^m J 

147 - 736m - 167m2 „ o\ , x (5267 + 10490m + 2387M2)7r2 



18mm ) ' ' 432m2 

291 - 959m (9 + 340m + 355m2 + ^2u^)t:'^ 56 



54m 18m^m 3 / 

(787 + 2607m + 1065m2 + 77M3)7r4 67^ 71387 , , , 
+ 720^15 + - T2% - 2A:2(«), 

, , , 4(1 - 3m + m2 + 5m3) , , 2, A 8(1 + 2m + 4m2)^ ^ ^ 
k2{u) = ^ ^^-^ (I2ni{u) + ln(2 - m) j + ^ ^S2,2(?^) 

2m -M^/ o 2t- / -^^ 2(1 + 10m + 4m2) 



3m3 



(2An2{u) - 27r2Li2(-M)) + J, ^ u{u) - 8 ln(M)Li3(M) 



5 - 10m + 16m2 - 4m^ . , ,2 8 + 36m - 33m2 - 20m3 . 

^ Li2(m)' + Li3(m) 

92 - 24m - 33m2 - 26m3 47 - 49m + 44m2 , , , 44, , 
+ Si,2(m) + ln(M)Li2(M) + — ln3(M) 

i^2 ll«„-,3 /I o«„-, A^,3\^2 



33 - 227m+ 157m2 - 116m^ (1- 26m -4m 



TT 



+ I 9^ + 3^1^ ) ^^^(^^ 

39 - 235m + 349m2 + 12MM7r2 (419 + 1178m + 815M2)7r2 

H ln2(M) - ^ ^ — 

18mm ^ ' 216m2 

807 - 2545m (12 + 65m + 41?/^ + 56^)712 A,,. 289 - 19m 



T^^^ + 14^3 ln(M) + Cs 

54m 18m^m / 18m 

m2 - 3 / . , , , ln2(M) + tt\ . \ 89437 
+ V '^"""^ ~ ln(M)Li2(-M) ln(l + m) j - 
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(227 - 1935m + 519^2 - 389#)7r^ {6 + 2u)tt^ 

1080^15 M2)-Co, 

, , , 10, , , 2 1 25 57r2 10^ 

, 8(6u + u^)^ 4(3 + Sm - 24^3 - 19^4) 144m + 28^3 
= 3m3 9^1^^ L^^(^) 9^1^ 

2(57 + 89m - 73^2 - 265^3 - 3mm V^) , , , 
+ 27^1^^ 

(23 + 315m - 507m2 + 4lM3)7r2 1111 - 6758m + 7951m2 
^ 54^13 + 162m2 ' 

, , , 4(13 + 8m)/ . , , , ^ -r. / N2 2t. / x\ 262 -298m, , , 
h{u) = -^—^^ ^ (2Lu{u) + 8S2,2(m) - Li2(M)2 - 7r'Li2(M) j + — ln(M) 

1-15m-67m2-3m3/ T- /-N 21 / n\ 4(1 + 21m+143m2 + 3iZ3)^ 
2Li3 M + tt' In m + ^ — ^Si,2 M 

4(2 - 12(7 - 29(72 - 3m=^) 5 - 171(7 + 507(72 + IGSm^^ 

ln(M)Ll2(M) Ll2(M) 

28m, , 12 - 41m + 190m2 + 91m^ , (3 + 214m + 35M2)7r2 

3(13 + 8M)7r^ 259 , , , 
+ 5m3 

fc6(«) = ^|^5^^^(l2^i(^) + ^'1^(2 - «)) - ^(24^2(^) - 27r^Li2(-M)) 
4(11 + 4m)_ ,_, 2(3 + 4m)/^ ^ . ,.2\ 13 - 56m - 8m2 



+ ^ ^3 ^ Li4(M) + (^8S2,2(«) - Li2(M)2 j ^^^35 In^(M) 

2(2 + 30m-38m2 + 3#),. 2(8 + 36m - 30m2 - IIm^) ^ 
H Li3(^i) H 7z Si,2(m) 

M'^M M-^M 

2(2 + 9u + 3m2) (6 + 53m + 81^2)772 269 
+ -I?. ln(ii)Li2(M) —^^ + 



M^M 



3u'^u 9m 



'2(11 -45m + 60m2 + 25m3) 2(17 + 12M)7r2\ ^ , , 2(7 + m), 

+ o... Li2(m) - ^-^Cs 



3m^m^ 3m3 

^2 _ Kn-,3\^2 



V? 



'203 - 257m (4- 18m-39m2 -5#)7r2\ , , , (155 + 108M)7r^ 
+ ( ^ TT^. — In M + ^ ' 



9m2 3m2m2 J ' ' 90m3 

2(8-6m + m2)/ 1 / NT- ^ N ln2(M)+7r2 



M^ 



(Li3(-M) - ln(M)Li2(-M) - ^"^^"^ ln(l + m)) , 



, , , 32^. 8(1 - 3m - 3m2 + m3)^ . 4(13 - 86m + 13m2) , , , 

fc,(M) = -Ll3(M) + ^-,3, Ll2(M) ln(M) 

^ 32(3 - M)7r2 _ 32^ _ 4(19 + 59m) 
9m^ m^ 9mm ' 
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ks{u) = ^ ^ (2Lu{u) + 8S2,2(m) - Li2(m)2 - 7r2Li2(M) j 

1 - 17m - 129^2 _ 101^3 _ 6^4 2i / x\ 28(1-2m) , 

+ (2Li3(n) + vr^ ln(n)) + ^3,, In^u) 

4(1 + 19m + 285^2 + 181^3 + 18^4) 3(19 + 32m + 12^2)71^ 

4(2 - 16m - 51^2 - 61m3) , , . (3 + 323m + 373^^ + 57u^)7r^ 
+ ln(«)Li2(M) + 

5- 187m + 363m2 + 1123^3 + 208m\ 262 - 771m + 149m2 ^ ^ 
+ 9^^ 

12 - 53m + 182m2 + 587m3 + 28m\ x 259 - 151m , , 

:r^3 1^ 5— 2A;9(m), 

3m"^M'^ 9mm 

fc9(M) = ^-^^ [I2niiu) + 7i' ln(2 - m) 

4(1 + 2m)/ ^, , , 2. / x\ 4(13 + 24m + 8m2)^ 
+ ^^^(247^2(u) - 2n^U2i-u)) - ^ ^-LUiu) 

2(5 + 8m + 8m2) / ^ ^ ^ T - . ^2^ 2(2 + 48m-28m2-25m3-6m4)^. , , 

- ^ 7 882,2 M - L12 M^ - ^ IT ^Ll3 M 

M* V / M*M 

2(8 + 70m + 16m2-71m^-14m^)„ , , 13 - 82m - 34^2 - SOm^ , 
2(2 + 17m + 17m2 + 6u^) (6 + 103m + 60m2 + 32u^)Tr^ 

F ln(M)Ll2(M) + 

M'^M 3M'^M 

/2(ll - 67m + 72m2 + 103m3 + 34m^) 2(23 + 40m + 24M2)7r2 ^ ^ ^ ^ 

_ /203 - 795m + 430^2 (4 - 26m - 47m2 - 103m3 - 2M^)7r2\ ^^^^^ 
V 9mm2 3u^v? / 

, 2(22 - 20m + 9m=^-2#) ^^. , wuw ^ ln\u) + 7r\ , , 

H 7 Li3(-M) - ln(M)Li2(-M) ln(l + 

\ 2 

2(23 - 17m + 10m2) (209 + 364m + 216M2)7r4 Stt^ 269 - 215m 

H 5 Cs T^f-^ in(2) —— , 

M'^ 90m* m 9mm 

, , , 32(1 + 2m) / . , , \ 8(1 - 5m - 21m2 - 17m3 - 2u^) ^ . 
hoiu) = ^(Li3(M) -Caj Ll2(M) 

4(13 - 124m - 223m2 - 38m3) , ^ ^ 16(15 + 4m + M2)7r2 
+ 9^^ 

^ 4(19 + 364m + 43m2) ^^^^ 
9m2m 

The definition of tlie functions 7ii,2(a^) and tlie constant Cq can be found in Section [3?T1 
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When we include a finite charm quark mass in our calculation, we encounter three 
additional coefficient functions which depend on the quark mass ratio z = rric/mf,. They 
have been introduced in P6l) and fHTl) and read 

^ /24m3(i - 2) - (49 + 119m + 175u^ - 55u^)z^ - 18(1 + u){l-8u + 3u'^)z^ 



+ 



9m2(1 + m)z2 

(69 - 80m + 27u'^)uz^ + 4(15 - 2u + 23u'^)z^ - 2(31 - \'iu)u^z'^ 



9m2(1 + m)22 



X ( Li,(.^) + 2 ln(.) 1„(1 4((3-19«y + 2(15-23S)S^^; . ^ 



+ 



9-u2'U 

4(24^3 - (69 - 80m + rju^)uz'^ + (49 + 119m + ITSm^ - SSm^)^^) 

9m2(1 +m)z 



/ . , ,,,,,, ,T , A 2(54M3-(21 -(37 + 37r2)M)Mz2 + 216Mz4), , , 
X ( Li2(-z) + ln(z)ln(l + z)j ^ riuuz'^ ^ ^^^^ 

_ 4:{l2u^ - (37 - 51m)mV + 2(15 - 32mm + u^)z^ - 9(1 - 8m + 3m2)(1 + u)z^) 

9m2(1 + m)z2 ^ 

2(72m3 - (469 + 37r2 - (59 - 37r'^)u)uz'^ + 12(29 + 3m + Uu^)z^) 4 

^^TTTT^V^ ^^^^^ ~ n^3 

27u{l + u)z^ 9 

(12(7 - m)m3 - lUu^z + 40(3 - m + 12^2)^4 _ 35(1 + ^)(i _ + 3^x2)^6)^2 

54m2(1 + m)2;2 

(6(69 - 80m + 27m2)mz3 - (163 - 189m)m222 _ 5(49 ^ ^^g^ ^ ^75^2 _ 55^3)^5)^2 

54m2(1 + m)2;2 

_ 108(27 + 11m)m3 - (5827 - 509m)m2^2 ^ 36(107 + 265m + 395m2 - 83u^)z'^ 

162m2(1 + m)22 

8((3- 19m)m2 + 2(15-23m)m22) 
+ TT^ Hz) ln(M), 

, , ^ 16^^^/ N 8(m2 + 4mz2) 8(m2-2m22)/ 

8(m3 - 2m2z + (I + 10m + m2);z2 _ 4(5 + 2m - m2)mz3 ^ ^(^i + m)(3 - m)mZ^) 

3m2m(1 + u) 

f T ■ / 2\ , , ^^ 2\\ 16((3 - m)m + (13 + 5m)mZ^) , , . 

X (L.,(.^) + 21n(.)Ml-.^)j 3„^(1+^) 
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— ^ ^ ^— ^- \Yy\z) ^ — ln(M) 

3u^u[\ + u) 9uu^ 

8(2^3 - 6u^z + (1 + 17^)^2 - 12(5 + 2u - v?)uz'^ + 3(1 + n)(3 - u)uz^)7i^ 

9u^u{l + u) 

16(^2 -2m22) A(U3 - 5u)u'^ + 6(86 + Alu - 13u^)uz^) 

ln[z) m[u) — 



3uu'^ 9u'^u{l + u) ' 

, , , 16(3^3 + (37+ 15m + 2m2)m22) 16(1 + 2m)^2 

kisiu, z) = ln{z) Ji(m, z) 

3'U^'u(l + u) 

8(1 -2n)(M2 + 4^^2) 16((1 -2m)u2 + 2(1 + 10u)u22) 



3uu{l + u) 3u'^u'^ 

32((2 -u + 3u^)u^z + (35 + 25m + 13^2 - u^)uz^) 
3m'^'u(1 + u) 

(t- f ^ , / M ^^ 8((1 -2m)m2-2(1 + 10m)m22w 
X { Li2(-z) + ln(z) ln(l + z)\ + ' ^ ^ f Li2(m) + ^'(m^ 

16(2(1 - 2u)u^ + (1 + 19m - 6u^)uz'^ + 3(1 + m)(5 + ^ 
^ 3u^u(l + u) ^ 

_ 8(1 - zf{{l - 2u)u^ - (5 - u)u'^uz + 3(1 + n)(5 + m)^^^) 

3m%(1 + m) 

4((43 - 29u)u^ + 6(140 + 127u + 79u^ - Au^)uz'^) 



+ 



(2(1 - 2u)u^ -3(2 -u + 3u^)u^z + (1 + 27m + 2u^ + 24u3)22)7r2 

9m%(1 + m) 

8(3(35 + 25m + 13^2 - m^)^^ _ 3(1 + m)(5 + M)2^)7r2 
9m3(1 + u) 

with the finite terms of the 4-topology Mis Xi 2(^5 z) that we introduced in Appendix [Cl 
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C Master Integrals 



In this appendix we present our results for the Mis that involve the charm quark mass. 
The Mis are normalized according to 



[dk] 



d/2 



with d = 4 — 2e. We suppress the H-ie-prescription of the particle propagators, 



Vi = {k- If - zV,, V^ = {p + kf - nf 

V2 = P- z^ml 



by 



V4 = {uq + ky 



and write z = rric/mb and u = 1 — u. 
3— topologies 

1 



[dk][dl] 




V1V2V3 



= (m 



2\l-2e 



J=-2 



with 






1 2 

= 2+^' 




5 , 
4 


Uq 


11 

= - + 6^ 



(3 - 2^2 + 2z^)7i^ 
6" 



(31) 



+ 6z^ + 2z^ (^(2 + z^) \n{z) - 7j \n{z) + 

{l-zY(u2iz^) + 2\n{z) \n{l-z^) 
55 15 



+ —z^ - — ( (8 + 12^2) \n\z) - 3(12 + 5z'^) \n{z) + (111 + in^) ] \n{z) 

X6 2 3 



(53) 



(54) 



(55) 



+ 2{l~z 



2'Hi{z) + Si,2(^') + ln(l - z^)U2{z^) + 2 ln(l - z)U2{-z) 

+ (2 ln2(l - z^) - ln2(l -z)- ln2(l + z)^ ln(z) - vr^ ln(l + z) + 2(3 

(15 + 24^ - 24^2 ^ 24^3 ^ 10^4)7^2 



+ 8^(1 + z") \hi2{-z) + ln(z) ln(l + z)J + 

_^,l_.,^(5.U..5.=,(L,(-).21nWMl-- 
and the function TCi{x) from Section [XT! 



12 



(56) 
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[dk][dl] 



[m 



2\-2e 



j:dr^e^ + Oie^)\ (57) 



J=-2 



with 








1 

~ 2' 






1 

~ 2' 




"0 




1 

2 " 


df'> 




11 




(3-^2)^2 

6 



11 0/ , , , \ , 2/ ^ (1 +4^ - 2^2)7^2 

- + 2^2 1^2 \n{z) - 3j 111^(2;) + ^ 



2ni{z) + Si,2(^') + ln(l - z')U2{z^) + 2 ln(l - z)U2{-z) 



+ %z\U2{-z) + ln(z) ln(l + 2) - (1 - 2)(1 + 3^) U2{z^) + 2 In(^) ln(l - z 



+ (^2 ln2(l - _ in2(i -z)- ln2(l + z)) In(^) - vr^ ln(l + z) + 2(3 
4— topologies 



(5^ 



The 4-topology Mis are complicated since they depend on three scales (Mm^,m^,m^). 
The divergent terms of the Mis are simple, 




[dk][dl] 



1 



V1V2V3V4 



[m 



2\-2e 



—, + [^- + lHu)\-^+X,{u,z) + 0{e] 



(59) 



and 



[dk][di\ 



V1V2V3V! 



(m 



^ - ( 2 In(^) + \n{u) ] — + l2{u, z) + 0{e) \ , 



u 



us 



(60) 
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but we could not find analytical expressions for the finite terms Xi 2('U, zj^. Unfortunately, 
a power expansion in z = mc/mf, cannot be applied for the given Mis, since it generates 
powers of z/{l — u) which invalidate the expansion for large values of u. 

We could, however, derive suitable one-dimensional Mellin-Barnes representations, 
which allow us to evaluate the finite terms numerically to very high precision. We use 
this representation for our numerical estimate of charm mass effects in Section |6l but also 
provide compact parameterizations which are in any case convenient since we expect the 
exact expressions to be rather involved. We propose the parameterizations 



Ii(u,z) 



I2(U,Z) 



Co — ClU 

1 + (1 -m)/'' 

^0 

{uo-uy{i-uy 



(61) 



which reproduce our numerical results over the entire u-range, < m < 1, and for physical 
values of 2; = mc/mf, at the percent level. The numerical values of the parameters 
(co, Cl, p, . . .) depend on the value of z and can be found in Table [TJ 



z 


Co 


Cl 


P 


do 


uo 


cr 


T 


0.250 


11.9022 


0.9512 


0.6575 


13.1958 


1.0164 


0.2077 


1.2244 


0.275 


11.5770 


0.9053 


0.6544 


12.3101 


1.0171 


0.2170 


1.2271 


0.300 


11.2538 


0.8523 


0.6506 


11.5442 


1.0177 


0.2256 


1.2295 


0.325 


10.9330 


0.7934 


0.6463 


10.8767 


1.0182 


0.2335 


1.2317 


0.350 


10.6147 


0.7293 


0.6414 


10.2786 


1.0187 


0.2412 


1.2339 



Table 1: Numerical input for the parameterizations ^61\) of the finite terms of 
the 4-topology Master Integrals. 
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